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The linear response of synchronized chaotic units with delayed couplings and feedback to small 
external perturbations is investigated in the context of communication with chaos synchronization. 
For iterated chaotic maps, the distribution of distances is calculated numerically and, for some 
special cases, analytically as well. Depending on model parameters, this distribution has power law 
tails leading to diverging moments of distances in the region of synchronization. The corresponding 
linear equations have multiplicative and additive noise due to perturbations and chaos. The response 
to small harmonic perturbations shows resonances related to coupling and feedback delay times. For 
perturbation from a binary message the bit error rate is calculated. The bit error rate is not related 
to the transverse Lyapunov exponents, and it can be reduced when additional noise is added to the 
transmitted signal. For some special cases, the bit error rate as a function of coupling strength has 
the structure of a devil's staircase, related to an iterated function system. Finally, the security of 
communication is discussed by comparing uni- and bi-directional couplings. 



I, INTRODUCTION 

Chaotic systems which are coupled to each other, can 
synchronize to a common chaotic trajectory [3, 0, E|. 
Even if dynamic systems exchange signals with a large 
delay time they still can synchronize completely, with- 
out any time shift 0, 0, [H, 0|. This phenomenon is 
presently attracting a lot of attention, partly because 
of its counter-intuitive fundamental aspect of nonlinear 
dynamics, and partly because of its potential for secure 
communication with chaotic signals [3, [3, [13], for a re- 
view see In fact, broadband communication with 
synchronized chaotic semiconductor lasers has recently 
been demonstrated over 120 km in a public fiber network 
EH- 

There are several ways to transmit a message with 
chaotic signals. One of them depends on a phenomenon 
which has been coined "chaos pass filter" [13, H] ■ It was 
observed that a chaotic system which is driven by the 
chaotic trajectory of its partner plus the (small) message 
responds essentially only with the trajectory but not with 
the message. Thus the chaotic system filters out any per- 
turbation and responds with its chaotic trajectory. Using 
this mechanism, the receiving unit can recover the mes- 
sage just by subtracting the incoming signal from its own 
variables. 

Adding a small message to a chaotic carrier has the 
potential for chaotic communication [llj. Obviously, it 
is difficult to extract the message from the irregular tra- 
jectory of the transmitted signal. However, when an ad- 
versary has access to the signal, and when he knows all 
the details including the parameters of the sending unit, 
he may be able to use an identical unit, sychronize it as 
well and extract the message. Recently, it has been sug- 
gested that even in this case secret communication may 
be possible. Two chaotic units which interact by bidirec- 
tional transmission have an advantage over an attacker 
driven by a unidirectional signal. This difference between 
bi- and unidirectional transmission opens applications on 



public cryptography fl^fTj^flT^. 

Chaos pass filter is observed in experiments on elec- 
tronic circuits and lasers and in simulations of chaotic 
systems. But it is not well understood. When chaotic 
units synchronize, their chaotic trajectories will be at- 
tracted by a synchronization manifold. Any perturbation 
of one of the two trajectories has two effects: Parallel to 
the synchronized manifold the perturbation will explode, 
since the system is chaotic. But perpendicular to the 
manifold the perturbation will decay to zero, thus the 
system will synchronize again. However, this does not 
necessarily mean that a permanent perturbation like a 
secret message is damped by the chaotic receiver. 

The linear response of chaotic systems to an exter- 
nal perturbation has been investigated before, mainly 
in the context of linear stochastic systems with multi- 
plicative and additive noise [13, [13]. In this paper the 
linear response of synchronized chaotic units is investi- 
gated. To obtain analytic results, we use iterated maps 
to study chaos pass filter, but many of the results are 
also observed in numerical simulations of chaotic differ- 
ential equations. In section UH the models are introduced 
and the linearized equations for the perturbations are de- 
rived. In section Hill we define quantities which measure 
the linear response to noisy signals. We observe large ex- 
cursions off the synchronization manifold which lead to 
a continuum of diverging moments. We calculate reso- 
nances, bit error rates, devil's stair cases and Lyapunov 
spectra. In section ITV[ a third unit, an attacker, is con- 
sidered which is driven by the exchanged signal and tries 
to synchronize with the communicating partners. The 
last section summarizes and discusses our results. 



II. THE MODELS 

Consider two chaotic units A and B which are cou- 
pled to each other by a function of their internal vari- 
ables. The exchanged signal arrives at the partner with 
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(a) Unit B is driven by unit A. 

a + m . 
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(b) Units A and B interact by mutual 
coupling. 

FIG. 1: Two different scenarios with two chaotic units. The 
states of units A and B are indicated by a and b, the exchanged 
small message by m. 



some delay time t. Since we want to compare our results 
with corresponding experiments on chaotic semiconduc- 
tor lasers, we add a self-feedback to each of the two units, 
again with a delay time 0, HO, HH • We compare two dif- 
ferent scenarios, indicated in Fig. [TJ (a) Unit B is driven 
by unit A, and (b) units A and B interact by a mutual 
coupling. For both cases, unit A is sending some message 
m to its partner B. This message is small and its content 
may be considered as small random noise. Thus we are 
interested in the linear response of unit B to noise. 

In the following, we investigate the two scenarios with 
chaotic iterated maps fillip. The dynamics is given by 

• case (a) of Fig. Q] 

a t = (1 - e)f(a t -i) + sf(a t - T ) 
b t = (1 - e)/(6 t _i) 

+ £Kf(b t - r ) + e(l - K)f{a t -T + mt-r), 



(1) 



case (b) of Fig. Q] 

a t = (1 - e)/(a t _i) 

+ £Kf(a t - T ) + e(l - K)f{b t - T ) 
h = (l-e)/(6t_i) 

+ SKf(b t - T ) + e(l - K)f(a t - T + m t _ r ). 



(2) 



There are two parameters: e measures the contribution 
of all delay terms while K measures the relative strength 
of the self- feedback. For the function / we mainly use 
two variants: 



• the Bernoulli shift 
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f(x) = —x mod 1 

and 

• the logistic map 

fix) = 4z(l - x) 



(3) 



(4) 



Note that all variables a t ,b t are iterated in the unit in- 
terval [0, 1]. 




FIG. 2: Phase diagram for Bernoulli units. Regions II + III 
(- -): synchronization of A and B for case (a), see Eq. (f5|); 
regions I + II ( — ): synchronization for case (b), see Eq. ©. 



Without noise, m t = 0, the spectra of Lyapunov expo- 
nents have been calculated for the Bernoulli system. In 
the limit of large delay, r — > oo, the phase diagram has 
been calculated analytically [12, US]- In the stationary 
state, the two units A and B are completely synchronized 
for 
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case (a), (5) 
case (b). (6) 



These regions are shown in Fig. [21 

The uni-directional system synchronizes in regions II 
and III and the bi-directional in regions I and II. Synchro- 
nization is complete, a t — bt, although the exchanged 
signal has a very long delay. 

The effect of small noise can be calculated by lineariz- 
ing equations |T]) and ([2]) in the vicinity of the synchro- 
nization manifold at = b t . The noise leads to a small 
deviation d t — b t — a t which is given by the linear equa- 
tions 



case (a) 
dt = (1 - e)fUdt-i 



case (b) 



+ £(1 - ^)lLr m t- 



(7) 



<k = (1 - e)f t _ 1 d t -x + £Kf't_ T dt-r 

+ e(l - «;)/t_ T • [mt-r - d t - T ) 
= (1 - e)fU<k-i + (2« - l)ef t _ T dt- r 

+ e(l - K )ft-T m t-r 



(«) 



where j[ is the derivative of f(x) at the synchronized 
trajectory a t — b t . Without noise, these equations de- 
termine the region of synchronization. For the Bernoulli 
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system one has f' t = 3/2 where we can ignore the jumps 
in the limit of infinitesimally small differences d t . In this 
case (and without noise) , the equations have been solved 
analytically and determine the phase diagram of Fig. [H 
In the synchronized regions, d t decays to zero whereas 
outside d t increases exponentially. 



III. LINEAR RESPONSE 

The transmitted message may be considered as noise 
with a low amplitude. Hence we want to understand how 
a synchronized chaotic system responds to small noise. In 
particular, we are interested in the linear response. 



A. Moments 

The distance from the synchronization manifold is de- 
scribed by the variable dt of equations (7|) and (JSj - We 
consider the regime of complete synchronization, i.e. pa- 
rameters for which, without noise, d t decays to zero. 
With noise, however, one obtains a distribution of dis- 
tances d t . We are interested in the stationary distribu- 
tion which develops after a transient time given by the 
transverse Lyapunov exponents. This distribution may 
be characterized by its moments. Hence, we define a re- 
sponse function 



Xn = , hm 



<M)^o (\m\ 



(9) 



where (. . .) means an average over the distribution of m 
and d. In the following we will use a symmetric distribu- 
tion of random numbers rtit , for example with a uniform 
distribution in the interval [—M,M]. 

Let us consider the simplest case, Fig. 1(a) with r = 1 
and k = 0. The chaotic unit A is driving unit B. The 
linearized equation (7J) takes the simple form 



d t = (1 - e)/^!^-! + ef t _ x m t -\ . 



(10) 



Let us assume that the variable f t -\ is an uncorrelated 
random number. Then this equation describes a linear 
process with multiplicative and additive noise. For the 
Bernoulli shift, f' t = § is constant and only additive noise 
remains. But for the logistic map, the distribution of /' 
is given by 



p(.n = 



i 



(11) 



It is well known that linear systems with multiplicative 
and additive noise may lead to stationary distributions 
which have a power law 0, [H, HH 



p(d) 



1 

dri 



(12) 



The exponent 7 is determined by 

(i-er- i <i/ ( 'r i ) = i 



(13) 



Hence, if the factor (1 — e)|/ t '| can take values larger than 
one, the equation lflO|) may drive the system to very large 
excursions, leading to a power law of the distribution. 
For the logistic map, the maximal slope is /' = ±4, thus 
we find large excursions of the variable 4 for e < |. 
We did not succeed to calculate the distribution of d t 
analytically, but the linear response can be derived. 

Squaring equation (fl0|l . one finds the second moment 
of the deviation from synchronization as 



(d 2 ) = j 



£ 2 (/ ,2 )(m 2 ) 

-(i-snn 



(14) 



For the Bernoulli shift, the second moment diverges for 
e < e c — I where synchronization does not occur. For 
the logistic map, however, the second moment diverges 
already inside of the synchronization region. Synchro- 
nization is given for 



ln(l 



<ln|/'|)<0 



(15) 



which gives e c — \. The second moment diverges 
when the denominator of equation (f!4|) becomes negative. 
Since one has (f 2 ) = 8, the second moment diverges for 



£ < £2 



-= = 0.646.. 

V8 



(16) 



The second moment diverges already inside the region 
of synchronization. That means that the distribution 
p(d) has developed a power law tail, and equation lflO|) 
describes rare large excursions of deviation d t from the 
synchronization manifold. Figure [3] shows a trajectory of 
dt for £ = 0.7 (b) where the second moment exists and 
for £ = 0.6 (a) where the susceptibility xi diverges. 

Power law tails of the distribution of a stochastic pro- 
cess with multiplicative and additive noise have been dis- 
cussed in the context of chaos synchronization, before 0|. 
This phenomenon has been called "on-off intermittency" 
and was discussed in the vicinity of the synchronization 
transition e c . Here the phenomenon is triggered by the 
message transmitted via the chaotic signal, and it is ob- 
served deep inside the region of synchronization. 

The linear equation iflOj) also determines higher mo- 
ments of p(d). Equation lfl3|l shows that Xn diverges at 
a coupling £„ given by the equation 



i=(i- e »ni/r> 



(17) 



For the Bernoulli shift, all moments diverge at the syn- 
chronization threshold e c , For the logistic map, however, 
all moments diverge at a different coupling £„ given by 



^r(f) 



(18) 
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FIG. 3: Deviation d t from the synchronization manifold for 
logistic map. (a) e = 0.6: large excursions, xa diverges; (b) 
£ = 0.7: X2 exists. Note the different scale in the second part 
of (b). 
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FIG. 4: The threshold e n below which Xn diverges. 



Decreasing e, the susceptibility Xoo diverges first. Then, 
at each value of e, a susceptibility with a lower value of n 
diverges until for n — > the synchronization threshold is 
obtained, where (ln|d|) diverges. Figure [4] shows e„ as a 
function of n. Note that Fig. |4] shows the exponent j(e) 
as well since one has 7 = n + 1. 

Figure [5] shows the power-law of the probability distri- 
bution p(d) of the distance d in a log-log plot for e = 0.65. 
In addition, the straight line corresponding to the expo- 
nent 7 is shown. It can be determined by Eq. (fl8|) (Fig. [4]) 
and 7 = n + 1 . 
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FIG. 5: Log-log plot of p(d), the histogram of the distances, 
for the logistic map, uni-directional scenario, r = 1, s = 0.65. 
Besides the small deviations from a perfect power-law for 
d « 0.01, there are two cut-offs: The first cut-off is for small 
d. It results from the additive noise which prevents that the 
distance very often becomes very small. The second cut-off 
is for large d. It results from the fact that \d\ is bounded 
between and 1 (since a and b are bounded between and 
1) which means that d cannot really diverge. 



U p to now we have discussed the driven system, 
Fig. l(a)| B ut the linearized equations for the interacting 
system, Fig. 1(b), are very similar: 



(1 - 2s)f t _ 1 dt-i + efUmt- 



(19) 



Hence, in this case, one finds diverging susceptibilities, as 
well. For the logistic map, the region of synchronization 
is given by 



1 3 

4 4 



and \n diverges at 
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?r(s) 



(20) 



(21) 



Fig. [6] shows X2 as a function of e. It diverges at the 
coupling strength e 2 = {\ ±575} = {0.323 . . . , 0.676 . . .}. 
With increasing coupling, the effect of the noise be- 
comes stronger, hence the susceptibility is not symmetric 
around e = 0.5. 

Up to now we have considered a system without de- 
lay. Now we are discussing the complete equations ([7]) 
and ([8]). If the delay time r is very large, the variables 
dt-i and d t - T may be considered as uncorrelated random 
variables. Hence the second moment diverges when the 
following inequalities hold: 



• Uni-directional 



l-[(l- £ ) 2 +eV](/' 2 ) <0 



(22) 
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FIG. 6: The susceptibility X2 as a function of £ in the case 
of uni- and bi-directional coupling and r = 1 for the logistic 
map. 




FIG. 7: Phase diagram for the unidirectional setup, the lo- 
gistic map and r — > oo, obtained from numerical simulations. 
The gray regime with black border ( — ) shows the synchro- 
nization region. The dashed line (- -) confines the region 
outside which X2 diverges. That is, the divergence already 
occurs inside the synchronization region. 



• Bi-directional 

l-[(l- £ ) 2 + £ 2 (2«-l) 2 ](/' 2 ) <0 



(23) 



For the Bernoulli shift one has /' = const. = § and 
all moments diverge at the synchronization threshold. 
For the logistic map, however, these equations cannot 
be solved analytically since the distribution of /' is not 
known. Figure [7] shows a numerical simulation which 
demonstrates that the susceptibility \2 already diverges 
inside of the synchronization region. 



B. Bit Error Rate 

In the following, we assume that A sends binary mes- 
sages (bits), i.e. m t = ±M with M <C 1 and (m) = 0. 
To reconstruct the small message m which A sends to 



B, unit B first has to subtract its own state from the 
received signal [ill ]: 



fh t = (a t + m t ) ~b t — m t - d t . 



(24) 



If units A and B were perfectly synchronized (d t = 0), 
the original message would be perfectly recovered by the 
above equation (m t — m t ). But since the original mes- 
sage is binary (±M), it is sufficient for a successful recon- 
struction if mt and rht have the same sign, m t m t > 0. To 
measure the quality of reconstruction, the bit error rate 
(BER) r is introduced, which is the ratio of the num- 
ber of incorrectly recovered bits to the total number of 
reveived bits: 



# incorrectly recovered bits 
# received bits 



(25) 



If the absolute value of the distance d is less than the 
absolute value of the messages, \dt\ < M, the message is 
recovered correctly with probability 1 because the sign of 
the message is not changed. In the other cases, \d t \ > M, 
the sign of the message is changed with probability i, i.e. 



(26) 



— M oo 
, =- I J p(d)dd + J p(d)dd 

V-oo M 

j P (d)dd\ . 



Note that this definition of the bit error rate assumes 
that the system has relaxed to a stationary distribution 
p(d). In fact, this definition may be only an upper bound 
since lower bit error rates may be achievable if one uses 
additional information about the transmitted signals [26| . 

The distribution of the distances, p(d), is known ana- 
lytically only in some special cases, see Appendix. In 
general, the bit error rate has to be determined by 
means of computer simulations. Figures El HO] and [TT] 
show the results of simulations for the uni-directional /bi- 
directional setup, for the Bernoulli map/logistic map and 
for no delay /large delay. 

For the Bernoulli map without delay, Fig. [HI we find 
a broad distribution of deviations d for small e in the 
case of uni-directional coupling, see Fig. fT27 a). With in- 
creasing e, Fig. [321(b), the values of the distribution get a 
fractal structure while the support stays connected. For 
e > | j Fig. [T2Tc). the distribution changes to a peaked 
structure with a fractal support. As derived in the Ap- 
pendix, this fractal structure leads to a devil's staircase 
for the bit error rate as a function of e. The bit error 
rate locks into rational values r = ^ with k and q being 



natural numbers and 



< r < 



In the bi-directional 



4 — ' — 2" 

case, the bit error rate is zero for | < e < | since for this 
interval the absolute values of the deviations are less than 
the message amplitude, \d\ < M, see Fig. [I27d). Simi- 
larly to the uni-directional case, one can find a staircase 
structure for the bit error rate for | < s < | . The fractal 
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FIG. 8: Bit error rate r for Bernoulli map, r = 1, k = 0, for 
uni- and bi-directional coupling. 




FIG. 9: Bit error rate r for logistic map, r = 1, K = 0, for 
uni- and bi-directional coupling. 



properties of the distribution p(d) are related to the the- 
ory of iterated function systems [27^. For the Bernoulli 
shift, Eqs. (fTO)) and (fl9|) give iterations of two linear func- 
tions with mt-i = ±1, Iterating a few randomly chosen 
functions can lead to fractal distribution. For details see 
Appendix. 

For the logistic map without delay, Fig. the distri- 
bution p(d) of deviations d has always a connected sup- 
port due to the multiplicative noise /', see Figures [T^e) 
and (f). The numerical simulations show a minimum of 
the bit error rate at e ~ 0.75 (uni-directional case) and 
e ~ 0.4 (bi-directional case). 

Note that the bit error rate is determined by the in- 
tegral Q26[l of the distribution of distances d from the 
synchronization manifold. We find that this integral is 
insensitive to long tails of this distribution, discussed in 
the previous subsection. Although the responses x« di- 
verge at different coupling strenghts e„, the bit error rate 
is smooth as a function of e. 

In the sequence of correctly and incorrectly recovered 
bits, there are correlations between time steps with a 



distance A of integer multiples of r and close to integer 
multiples, i.e. A w kr with k = . . . , —2, —1, 0, 1,2,.. ., 

e.g. A = 1, A = 2r — 2, An example can be seen in 

Fig. [13l For short time shifts A, the bit errors are not 
correlated. But this does not necessarily mean that the 
bit error rate can always be reduced by using larger bits 
(bits longer than one time step), since the system imme- 
diately responds to a long constant perturbation, i.e. the 
errors each time step are not independent anymore if the 
bits are longer than one time step. However, in most of 
the cases the bit error rate is reduced by using longer 
bits, as seen in Fig. [T4l 



C. Noise-reduced bit error rate 

In the previous subsections we have investigated the 
case where a message m t — ±M is added to the signal 
which unit A sends to B, and we have calculated the 
response of unit B. This response was quantified with 
the bit error rate r B . 

Now we extend this scenario to the case where unit B 
adds noise w (or a message) to its transmitted signal as 
well, 



a t = (1- e)/(a t _i) 

+ £nf(a t - T ) + e(l - K)f(b t - T + w t - T ) 
6 t = (l-e)/(fc t _i) 

+ enfibt-r) + e(l - «)/(a t _ T + m t _ T ) 



(27) 



see also Fig.HU Surprisingly, we find that this additional 
noise can reduce the bit error rate r B . 

This effect can be seen with the Bernoulli map as well 
as with the logistic map for many of the possible choices 
of the parameters in Eq. (|27|l . However, to treat the 
effect analytically, we set r = 1, n = 0, e = |. The 
corresponding equations are 



a t = 2/K-1) + \^ ht - x + 



(28) 



Let us consider the case of the logistic map. Obviously, 
for the above choice of parameters A and B synchronize 
immediately when the message m t and the noise w t are 
switched off. When message and noise are small, the de- 
viation from the synchronous trajectory is obtained from 
the linearized equations. With b t = a t + d t one finds 



7 ft—i 1 
d t = ——{m t - 



Wt-l) 



(29) 



where f' t _ 1 is the derivative of the nonlinear function / 
at the perturbed trajectory at; it may be considered as a 
random number. In the case of the logistic map, Eq. ((H), 
the distribution of s = f'/2 is given by 



fa) 



1 



7iV4 - s 2 



(30) 
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FIG. 10: Bit error rate r (gray-scale) for Bernoulli map, r = 50. The dashed lines (- -) represent the borders of the 
synchronization regions, see Fig. [2] The white dotted line (■ ■) in (b) confines the region where r < 10 -4 . 
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FIG. 11: Bit error rate r (gray-scale) for logistic map, r = 50. The bands around £ ~ 0.15, 0.45 and 0.55 represent periodic 
windows. 



The unit B recovers the message from the sign of 



The bit error rate of B is given by 

?* B = prob(r7Jt mf < 0) 

= prob(Af 2 - St-iM 2 + st-iWt-WH < 0) 
= prob(M - st-iM ± st-iivt-i < 0) 
= prob(M - s(M + w) < 0) 



Pw(w) / p(s) ds dw 



Pw{w) 



— oo 
-|M 



p(s) ds dw 



-see -1 ! 2 



w 

1+ M 



m f = a t + m t - b t — m t - d t 



(31) 



p w (w) dw 
(32) 
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FIG. 12: Probability distribution (in arb. units) for the deviations d with r = 1, k = 0. 
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FIG. 13: Auto-correlation C with time shift A of the sequence 
of correctly and incorrectly recovered bits, C = ^ggyzx^ i 
c t = +1: correctly recovered bit, c< = —1: incorrectly re- 
covered bit. Logistic map, bi-directional scenario, r = 50, 
£ = 0.7, K = 0.6. 



For small values of the noise, jj <C 1, one can expand 
the inverse sec; for a symmetric distribution of noise one 
finds 



r B = i _ yy) 

3 6V3M 2 tt 



205(w 4 



108\/3M 4 7r 



0((w 6 )) (33) 



For large noise we have calculated the integral l(32|) for 
a flat distribution, w S [— W, W] with p w (w) — . The 
result is given in Appendix [Bj and Fig. [16] shows the 
bit error rate r B (W) as a function of noise W. If unit 



FIG. 14: Bit error rate (•) as a function of the bit length for 
logistic map, bi-directional scenario, r = 50, e = 0.7, k = 0.6. 
For bit lengths > 1, the bit error rate deviates from a behavior 
(□) of independent bit errors (binomial sum). 




a + m 




B 



b + w 

FIG. 15: B sends noise w in addition to the signal b. 



B adds noise to its transmitted signal, it reduces its bit 
error rate. Only when the ratio between maximal noise 
W to the strength of the bit message M is larger than 
the value « 3.4, B has a larger error rate than without 
noise. 
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0.5 r 



0.4- 




0.2 



0.1 



W/M 



FIG. 16: Bit error rate r B for different strengths W of the 
noise w £ [—W,W] with flat distribution p m (w) = ^7 which 
B sends in addition to the signal b. Surprisingly, the bit error 
rate decreases due to the noise. The analytical result for the 
logistic map with t = 1, k — and e = \ is shown. The 
horizontal dashed line represents the bit error rate without 
noise, W = 0. 



W/M 

FIG. 18: Bit error rate r B for different strengths W of the 
noise w € [— W, W] with flat distribution p w (w) = ^4?. Result 
from simulation for logistic map, e = 0.9, K = 0.4, r = 50. 
The dashed line represents the bit error rate without noise, 
W = 0. 



0.5 r 




The main idea of this explanation for the special case 
also holds for the general case of Eq. l(27|): If the prob- 
ability density p(d) without noise prevails in the region 
\d\ > M, changing d by noise leads to an increasing of 
the probability density within the interval \d\ < M which 
decreases the bit error rate. 



0.2 



D. Resonances 



0.1- 



W/M 



FIG. 17: Bit error rate r for different strengths W of the 

2^7. Result 
ax mod 1 with 



noise w £ [— W, W] with flat distribution p w (w 
from simulation for Bernoulli map f(x 

a = § , e — 0.95, K = 0.45, r = 50. The dashed line represents 
the bit error rate without noise, W = 0. 



As mentioned before, the effect of decreasing the bit 
error rate by noise exists not only for the special case of 
Eq. i|28j) but also for the more general case of Eq. f27j) 
with delay and feedback (in certain regions of the pa- 
rameter space), both for the Bernoulli map and for the 
logistic map. This is shown in Figures [T7l and [T8l 

In order to explain the positive effect of noise, let us go 
back to the special case of Eq. I|28p and the logistic map. 
Consider the case where there is a bit error without noise. 
According to Eq. (J29j) , with w t = 0, one has \f'\ > 2. In 
this case noise can reduce the value of d t and remove 
the bit error. Noise can also generate a bit error when 
I /'| < 2. But since the distribution of |/'| is increasing, 
the probability of the first case is larger, resulting in a 
reduction of the total bit error rate. 



In the previous subsection we have investigated the 
linear response of synchronized chaotic units to a ran- 
dom perturbation m t . The response was described by the 
functions Xn, Eq. (J9j) . For the Bernoulli map, the corre- 
sponding equations lO and iJH) are linear equations with 
constant coefficients. Hence, in this special case, the re- 
sponse of the chaotic system can be described in terms of 
harmonic perturbations m t — Mexp(— iu)t). The system 
responds with the identical frequency d t = D exp(-iujt), 
and the complex amplitude D is given by 



uni-directional case 



D = M 



e(l — n)a 



exp(iu>T) — (1 — e)a exp(zw(r — 1)) — sua 



(34) 



• bi-directional case 



D = M 



s(l — K)a 



exp(io;T) — (1 — e)a exp(icj(r — 1)) 



e(2k- l)a 
(35) 

with a = f = |. The unit A sends the harmonic signal 
mt and the unit B tries to recover it by subtracting its 
own variable bt from the received signal a t + m t . The 
recovered signal is m = at + mt — bt = mt — dt- Hence 
the amplitude of the recovered signal is given by 



\m\ = \M - D\ 



(36) 
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FIG. 19: Amplitude of the reconstructed signal for different 
signal frequencies. Bernoulli map, uni-directional coupling, 
e = 0.8, k = 0.55, r = 50. 




2/10 3/10 
frequency v 



5/10 



FIG. 20: Fourier transformation of the reconstructed signal, 
S( v ) = E^Li^texpt^t- l)(Ni> - 1)]. Logistic map, uni- 
directional coupling, e = 0.7, k = 0.4, r = 50, lo — 2tv ■ 
-. In the Fourier spectrum there is a strong peak at v = 
Yq = - which corresponds to the frequency of the harmonic 
perturbation. 



Figure [19] shows |m| as a function of u> close to the 
phase boundary of synchronization. The amplitude of 
the restructed signal shows peaks (resonances) and/or 
cuts with distance — . The resonances diverge at the 
synchronization boundary. 

For the logistic map the coefficients of the linear equa- 
tions (7]) and ([8]) depend on time. Thus an exact Fourier 
decomposition of the transmitted signal is not possi- 
ble. However, our numerical results of Fig. l20l show that 
the corresponding Fourier component of the response rht 
shows resonances as well [28]. 

These results show that a chaotic system can func- 
tion as a sharp harmonic filter. The transmitted signal 
is irregular, and the harmonic perturbation is arbitrar- 
ily small. Nevertheless, the reveiver can filter out this 
perturbation with high precision. It would be interesting 
to investigate this harmonic filter with respect to secret 
communication. In fact, the response of synchronized 
chaotic semiconductor lasers to a harmonic perturbation 



has been investigated in Ref. [2£ 



E. Transverse Lyapunov Spectra 

In the previous subsections we have investigated the 
response of a synchronized chaotic system to persistent 
random or harmonic perturbations. The chaotic dynamic 
tries to bring the system to the synchronization manifold, 
whereas the perturbation drives the system away from 
this manifold. The competition between these two mech- 
anisms results in the linear response described above. 

The first mechanism - relaxation to the synchroniza- 
tion manifold - is usually described in terms of transverse 
or conditional Lyapunov exponents. In the synchronized 
region those exponents are negative and the largest one 
becomes positive just at the phase boundary, resulting in 
a slowing down of the relaxation to synchronization. As 
we have seen above, at the phase boundary the response 
to perturbation diverges, as well. 

Following these arguments, one might expect that the 
qualitative behavior of the response functions %„ or the 
bit error rate r is related to the behavior of the largest 
transverse Lyapunov exponent A tran s- However, this is 
not true. For example, consider the simple driven sys- 
tem without delay, Eq. (fTO]) , The transverse Lyapunov 
exponent is given by l(l5|). 



ln(l 



(37) 



Thus, it decreases monotonically with increasing coupling 
strength e. However, the bit error rate r first decreases 
with £ to a minimum value but then it increases again. 
This holds for the Bernoulli as well as for the logistic 
map, see Figs. [8] and [9] 

With delay r, the spectra of Lyapunov exponents can 
be calculated by solving polynomial equations of order r 
[13]. In Fig. [21] the spectrum of r transverse exponents 
is shown as a function of e for the self-feedback strength 
n = 0.2. Again, the qualitative behavior of the largest 
exponent is not related to the one of the bit error rate r. 



IV. ATTACKING WITH DRIVEN UNITS 

The main motivation of this investigation is the appli- 
cation of chaos synchronization to secure communication. 
Let us assume that two partners A and B have realized 
a synchronized chaotic system and transmit messages, as 
described above. But now an attacker E has access to the 
transmitted signal. We assume that the attacker can only 
record the transmitted signal, but he cannot send his own 
signal to the partners A and B. Thus E is driven by the 
transmitted signals but he cannot interact. This impor- 
tant difference between partners and attacking units has 
been suggested for a novel public key encryption protocol 

am- 

To extract the message from the transmitted signals, 
the attacker has to synchronize with the two partners. 
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FIG. 21: Spectrum of transverse Lyapunov exponents A and 
bit error rate r for Bernoulli map, uni-directional coupling, 
k — 0.2, t = 50. Except for the synchronization bound- 
ary where the largest transverse Lyapunov exponent becomes 
negative and the bit error rate becomes less than i the qual- 
itative behavior of the Lyapunov exponents and the bit error 
rate are not related. 



Thus it is obvious that he will use the same chaotic units 
with identical parameters. If he does not know those 
parameters and if he is not able to estimate them from 
the data, he will not be able to synchronize. But se- 
cure cryptography does not use any secret parameters, 
hence we assume that the attacker E knows everything 
wha t B k nows about A. In this case, the configuration of 
Fig. 1(a) A is driving B, is obviously useless, since E is 
in the same situation as B and can restore the message 
of A with identical precision a s B. 

Thus we consider only Fig. 1(b) where A and B are 
mutually interacting. The attacking unit E may have 
access to both transmission lines, as indicated in Fig.l22l 
The corresponding equation for the units A and B are 
given by Eq. |(2|), and E uses the equations: 



e* = {l-e)f{e t -i) + e<7f{b t -r) + e{l-(j)f{at-T + mt-r) 

(38) 

If unit E adjusts its parameter to a = k, it is obvious that 
unit E receives the identical drive including the same self- 
feedback as unit B . Consequently, E achieves an identical 
bit error rate. 

The situation is different if B sends a message or noise 
as well. This scenario was investigated in Section UlI CI 
In this case and for t — 1, the equations for A, B and a 
simple attacker E are the following: 

at = (1 - e)/(at_i) + ef(b t -i + w t -i) 

b t = (l- e)/(6t_i) + E/(a*_i + rm-i) (39) 

e t = (1 - e)f{b t -x + w t -i) + ef(a t -i + m t -i) 

From the above equations, the distance between A and 
B, dt, and the distance between A and E, gt, can be 



calculated: 
d t = b t ~ a t 

= (1 - 2e)f f _ 1 d t - 1 + e/j^mt-i - ef^wt-i 
g t = e t - a t 

= (1 - 2e)fUdt-i + efUmt-i + (1 - 2e)f t _ lWt -i 

(40) 

If one sets e = \, the bit error rate of B is reduced by 
the noise w, see Sec. MI CI whereas the bit error rate of 
E is then independent of w and is not reduced (indicated 
by the dashed line in Fig. [T6|) , This phenomenon has 
been reported for chaotic semiconductor lasers. A large 
difference in bit error rates for the partners and for the 
attacker was observed for the Lang-Kobayashi equations 
describing chaotic semiconductor lasers. The effect was 
named "mutual chaos pass filter" [l(| • 

However, for our simple model, an advanced attacker 
can achieve the same bit error rate as B. This should 
be shown for arbitrary e. The equation of the advanced 
attacker is: 



e t =(1 - 2e)/(e t _i) + ef(a t -i + m t - 
+ £f{b t -i + u>t-i) - efi-iVt-i 



l) 



(41) 



Here, y is an additional noise with the same distribution 
as w, and f' t _ 1 can be calculated from cit-i and bt-i. 
For other maps than Bernoulli map and logistic map this 
might not be possible. Then it is sufficient to use val- 
ues with the same distribution as /'. For the advanced 
attacker the distance between A and E is: 



9t 



(1 - 2e)fl_ l g t - 1 + sfUm-i ~ efl-iVt-i (42) 



which is identical to the equation for d t and thus results 
in the same bit error rate. 

An attacker is also successful if r > 1. The dynamics 
of units A and B and of an advanced attacker E are: 

a t = (1 - e)f(a t -i) + enf{a t - T ) 

+ £(1 - K)f(b t -T + W t -r) 

bt = (1 - e)/(&t_i) + enfQh-r) 

+ e(l - K)f(a t - T + m t _ T ) 
et = (1 - e)f(et-i) + s(2k - l)/(e t _ T ) 

+ e(l - K)f(at-T + mt-r) + e(l - «)/(&t-r + m-r) 

- e(l - K)ft- T Vt-T 

(43) 

So the distance of the attacker to unit A, gt = e t — a t , 
corresponds to the difference between A and B, dt = 
b t - a t : 



dt = (1 - e)ft_ x dt-i + e(2n - l)ft_ T d t - r 
+ e(l - K)f't_ r {mt- T - w t -r) 

9t = (l- e)fl-i9t-i + e{2n - l)#_ rSt _ T 
+ e(l - K)/t-T TO t-r - e(l - n)f' t _ T yt-T 



(44) 



Again, if the attacker E is able to add additional noise y t 
where f' t yt has the same or similar distribution as f' t wt, 
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a+m 



a+m 




B 



E 



FIG. 22: The two units A and B are coupled as in Fig. 1(b) 
Here, an attacker E has access to both signals. 




A 



a+m 



a+m 




B 



FIG. 23: The two units A and B are coupled as in Fig. 1(b) 
Here, an attacker E has access to the signal which goes from 
A to B. 



he can reduce his bit error rate in a similar way as partner 
B. 

Note that our results are based on iterated maps. For 
more complex systems, like the Lang-Kobayashi equa- 
tions and other higher-dimensional chaotic systems, it is 
not obvious that the same conclusions about the mutual 
chaos pass filter hold. 

Now let us come back to the general case r > 1 and 
let us assume that the attacker E is only able to record 
the transmission from A to B, see Fig. [23l For this case, 
the equation for E is 

e t = (l-e)/(ei_i)+ecr/(et-r)+e(l-o-)/(at-T + mt-r) 

(45) 

We want to compare the linear response of unit E to the 
one of B. According to the phase diagram Fig. [21 without 
message, unit E can always synchronize to A and B if 
it uses an identical parameter e and keeps its feedback 
parameter a in region II or III. 

Let us assume that A and B take their parameters 
(e, k) in region II. If the attacker E chooses k = <r, he 
obtains the same response and bit error rate as B. 

Figure [24Ta) shows the second moment and bit error 
rate of E as a function of a. Both quantities decrease 
monotonically with the feedback parameter ex. This re- 
sult may be unexpected. One might expect that the at- 
tacker E obtains the lowest bit error rate if he uses a 
strong signal containing the message. The opposite is 
true, E can even achieve a lower bit error rate than B if he 
takes the almost weakest possible contribution of the sig- 
nal without losing synchronization, here a « ^ = 0.583 
for e = jjy — 0.8 according to Eq. ((!]). Only very close to 
the synchronization boundary the response xf diverges 
and the bit error rate increases to its maximal value 
r E = 0.5. 

If the two partners use parameters of region I, the at- 




0.5 
0.4 
2 0.3 
0.2 
0.1 




FIG. 24: Bit error rate r E and response function xf as a 
function of a, for Bernoulli map, r = 50, e = 0.8 and Eqs. (fSJ) 
and l)45p. The values for unit B (which are not dependent of 
a) are shown as dashed lines, (a) k — 0.4 ; (b) « = 0.75 



tacker E can never achieve the same response as partner 
B. Figure l24Tb) shows that the bit error rate as well as 
%2 of E is much larger than the ones of partner B. The 
bit error rate of E is always at least about 10 times larger 
than the one of B. 



V. SUMMARY 

Chaos synchronization has been investigated in the 
context of secure communication via chaotic signals. In 
particular, when a small secret message is a perturbation 
of the sending chaotic unit, the receiver can recover this 
message by a mechanism which has been named chaos 
pass filter. Thus chaos pass filter is the linear response 
of a synchronized chaotic system to a quasi-random bi- 
nary perturbation. This perturbation has a component 
which is acting transverse to the synchronization mani- 
fold. Hence the message drives the system away from syn- 
chronization whereas the dynamics of the system, quan- 
tified by conditional Lyapunov exponents, brings the tra- 
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jectories back to the manifold. The competition between 
these two effects leads to the chaos pass filter which allows 
the receiving unit to recover the message by subtracting 
its own chaotic trajectory from the received signal. 

Consequently, one expects that perturbations of the 
sender are just damped by the receiver; the receiver fil- 
ters out the perturbation and responds essentially with 
the unperturbed trajectory. However, our numerical and 
analytical investigations of iterated maps show that the 
mechanism of chaos pass filter is much more complex. 
The response of the receiver to the perturbation of the 
sender can be very large. Close to the synchronization 
boundary it diverges. But even deep inside the region 
of synchronization where the transverse Lyapunov expo- 
nents are negative and large, the receiver makes large 
excursions away from the synchronization manifold re- 
sulting in a power law and diverging moments of the dis- 
tribution of deviations between the sending and receiv- 
ing trajectory. Mathematically, this behavior is a conse- 
quence of multiplicative and additive noise appearing in 
the equations of linear response. 

In addition, the linear response to a periodic pertur- 
bation shows resonances due to the delayed self-feedback 
of the sending unit. Depending on parameters and fre- 
quency, those resonances can be very large but the re- 
sponse can also be suppressed. 

A message produces a broad distribution of deviations 
between the trajectories of sender and receiver. The bit 
error rate is given by an integral over this distribution. It 
turns out that the bit error rate approaches continuously 
the maximal value 50% at the synchronization boundary, 
as expected. But inside the region of synchronization the 
bit error rate shows a complex nonmonotonic behavior 
which cannot be related to the properties of transverse 
Lyapunov exponents. 

In special cases we could calculate the bit error rate 
analytically. Relating it to an iterated function system 
we found a fractal distribution of deviations yielding a 
devil's staircase for the bit error rate as a function of 
model parameters. 

Adding additional noise to the chaotic transmitted sig- 
nals can reduce the bit error rate. This opens the pos- 
sibility to improve the security of chaos communication 
with bi-directional transmission of chaotic signals. When 
the receiver adds an additional message or noise to its 
transmitted signal he can improve his bit error rate (mu- 
tual chaos pass filter). 

With respect to secure communication, the difference 
between uni- and bi-directional couplings of the two part- 
ners has been investigated in detail. Assume that a third 
attacking unit has knowledge about all details of the 
chaotic systems and can record the transmitted signals. 
When the two partners are coupled uni-directionally, a 
third attacking unit can synchronize as well, and retrieve 
the secret message. However, when the two units are in- 
teracting by a bi-directional transmission of signals, an 
attacker may have problems to synchronize as well. It 
turns out that an attacker who can record and adjust 



both lines of transmission can synchronize. In this case 
even the mutual chaos pass filter does not help and one 
needs more advanced methods like secret commutative 
filters to establish a secure communication [17]. But 
when the attacker can only adjust to one direction of 
transmission, one can find parameters such that the bit 
error rate of the attacker is much larger than the one of 
the partners. 
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APPENDIX A: ANALYTIC RESULTS FOR THE 
BIT ERROR RATE 

Outside the synchronization region the bit error rate 
trivially equals h. Inside the synchronization region in 
general, the bit error rate has to be determined by means 
of computer simulations. However, for some special cases 
it can be calculated analytically. 



1. Logistic map, uni-directional setup, r = 1, e — 1, 

K = 



In the case of the uni-directional setup (case (a) of 
Fig. [TJ with the logistic map and no delay (t = 1), the 
bit error rate can be calculated for the point e = 1. The 
dynamics is given by 



f{fl t -x) 
f(a t -i + m t - 



from which follows that 



dt = ft-imt-i 



(Al) 



(A2) 



From the facts that (1) f' t _ l and m t -i are uncorrelated, 
(2) the probability distribution of /' is symmetric about 
/' = and (3) mt — ±M, follows that d has basically 
the same probability distribution as /'. Therefore, the 
bit error rate can easily be calculated, see Eqs. f26|) and 
CD: 



arcsec(4) 



0.4196 . 



This is in agreement with the numerical simulations 
shown in Fig. [9] 
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2. Logistic map, bi-directional setup, r = 1, e = 5, 

K = 



Similarly, in the bi-directional case, the bit error rate 
can be calculated for e — \: The dynamics is given by 



a t = 7;f{a t -i) + ^f(h-i) 

b t = -^fih-i) + g/K-i + m t-i) 



(A4) 



from which follows that 



d t = o/t'-l" 1 *-! 



Comparing this with Eq. (| A2|) leads to 

P(j')d/' 




(A5) 



(A6) 



which is also in agreement with the numerical simulations 
shown in Fig. [9] 



indicated by gray stripes). This fact makes the distri- 
bution p(d) complicated. How ever, f or e > | the two 



maps have no overlap, see Fig. 26(c) and the distribu- 
tion is manageable analytically. As a result of the gap 
in the co-domain (indicated by a zigzag pattern), a gap 
in the domain emerges in the next time step. The latter 
gap produces two further gaps in the co-domain which 
become gaps in the domain in the next time step. The 
result of this iterative process is that the distribution 
p(d) has a fractal support. The first and largest gaps are 
shown in Fig. [27] The first gap is called G. The gaps 
produced by G are G~ and G + . The gaps coming from 
G + are called G h and G ++ ; the gaps coming from G~ 
are called G and G H and so on. 

Now we want to calculate the exact position of the 
gaps. The fixed points of dr and d + are called d~ and 
d+ . One can easily calculate that 



dt 



3eM 
3e-l 

3eM 
3e- 1 



and 



(A8) 
(A9) 



From Fig. 26(c) it can be seen that 



3. Bernoulli map, r = 1 

In parts of Fig. [8] one can discover a staircase structure 
for the bit error rate. For the uni-directional setup this 
is true for e > | , while for the bi-directional setup this is 
valid for | < e < |. For these regions the bit error rate 
can be calculated analytically. If one takes a closer look 
at the staircases, Fig. [25] it becomes apparent that they 
have infinitely many steps, i.e. they are a kind of devil's 
staircase. 



G 



d-(dt),d+{d-)[ 

3e(3e-2)M 3e(3e-2)M 



3e- 1 



3e - 1 



(A10) 



which is about ]-0.19M, 0.19M[ for e = 0.7, see Fig. [27] 
The gap G + is generated by applying d + to G, i.e. 



G+ =}d+(d-(df)),d+(d+(d-))[ 

3e(5 - 12e + 9e 2 )M 3e(5 



12e + 9e 2 )M 



2(3e - 1) 



2(3e- 1) 



(AH) 



a. Uni-directional coupling 

This staircase structure should be explained here for 
the case of uni-directional coupling. The equation for 
the distance can be written in the following way, see also 
Eq. CED: 



dt = 



d t = -(l-e)d t - 1 
# = f(l-e)dt-i 



(A7) 



The two equations represent the two different bits. If d t 



is plotted versus dt— 1, then d t and 



d( are two paral- 



lel straight lines, see Fig. [261 The values of this itera- 
tion dt{dt-i) generate the distribution p(d) from which, 
in principle, the bit error rate can be calculated. The 
dashed boxes ( — ) in Fig. [26] represent the interval which 
d t is bounded to (due to the attracting fixed points). For 
e < I the two maps d^ and df ha ve a certain overlap 
in their co-domain, see Fig. 26(a) (the co-domains are 



which is about ]0.96M, 1.14M[ for e = 0.7, see Fig. EH 

Then the gap G h , for example, is generated by ap- 
plying d~ to G + , i.e. 

G~ + =]d-(d + (d-(df))),d-(d+(d + (d-)))[ (A12) 

and so on. 

Due to the constant and equal slope of d~ and d + , 
and due to the fact that there is no overlap between the 
co-domains of d~ and d + , the relative frequency of all 
distances d which occur (i.e. which are not inside a gap) 
are equal. This means that Fig.[27|can also be seen as the 
corresponding histogram; all bars have the same height. 

The bit error rate is related to the integral from — M to 
+M over the distribution of the distances, see Eq. (f26|) . 
From Fig. [27] it can be seen that for e = 0.7 this integral 
exactly equals \; thus, the bit error rate equals \, which 
is in agreement with Fig. [25] 

If e is changed, then the positions of the gaps are 
changed, too. As long as the gap G + contains the value 
+M (= as long as the gap G - contains the value — M), 
the integral yields \ and the bit error rate is \. This 
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PIG. 25: Bit error rate r for Bernoulli map, r = 1, K = 0, uni-directional coupling, see also Fig. Zooming in the staircase 
structure reveals more and more steps. 




(a) £ = 0.6; d and d + have an overlap (b) s = | (c) £ = 0.7; d and d + have no overlap 

in their co-domains. in their co-domains. 



FIG. 26: The iteration d t (d t -i) = df(d t ~i) of the distances for different coupling parameters e. The dashed boxes represent 
the interval which d is bounded to due to the fixed points. Additionally, the bisecting line d t (d t -\) = dt-i is plotted. The gray 
stripes show the co-domains of d~ and d + . 



explains the plateau AB in Fig. EH With the aid of 
Eq. IjAlip we can calculate the exact position of this 
plateau: 

2 

A: e = - = 0.6, (A13) 

B: e= -(1 +v / 2) » 0.804738. (A14) 
Similarly, one gets the point C of Fig.[25l The bit error 



rate becomes \ when the integral starts to be 0. This is 
when the gap G is as large as (or larger than) the interval 
[-M,M]. Considering Eq. ([ATO]) yields: 

C: e = -(3 + V5) » 0.872678. (A15) 
6 

The plateau DE, which has the value I, can be calcu- 
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k=Q 



k= 1 



k = 2 



k=3 



k=4 




-2M -M 







M 2M 



• f < W < |M 



n %= I sec" 1 (2 + dw (B2) 



• W >\M 



2irW 



W r W-2M 



M/2 J-M/2 



sec 1 2 + — dw 



The result is: 



• W < 



M 



r B = 9(W) - g(-W) 



(B3) 



(B4) 



FIG. 27: Gaps in the domain/co-domain/distribution of d 
and d + for different recursion depths k. s = 0.7 . 



lated considering the gap G + . One gets: 

D: £ w 0.837266 (A16) 
E: £ re 0.866386 (A17) 

All other plateaus can be calculated with the aid of 
smaller gaps. 



b. Bi-directional coupling 

The calculations for the case of bi-directional coupling 
are very similar to the ones for the uni-directional case. 
Here, only few results should be shown. 

For 



- 1 5 - 

0.3 = - < e < - = 0.5 

3 ~ ~ 9 



the bit error rate is 0. 
For 



11 2 
0.5749 «- + ^<£<-=0.6 
6 V6 ~ ~ 3 



(A18) 



(A19) 



• f < W < |M 



• W > IM 



r B = g{W) 



r B = g(W) + g(W - 2M) 



with 



(B5) 



(B6) 



2 



„ 2w 

2 H 

M 



„ 2w , 

2 H 1 

Af 1 



(B7) 



the bit error rate is \. 



APPENDIX B: NOISE-REDUCED BIT ERROR 
RATE 

The integral of Eq. j[32|) can be calculated if the noise 
has a flat distribution in the interval w G [-W, W] with 
p w (w) = • The problem can be split up into three 
intervals. 

• W < § 

r B = l W sec" 1 (2+^] dw (Bl) 

2nW J_ w \ MJ y ' 
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